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Abstract 

Assuming the idea of retardation as an underlying axiom, we investigate how Jefimenko's and 
MaxweU's equations can be inferred. In the inference, we begin with the retarded versions of 
Coulomb's and Biot-Savart's field expression as an incomplete, starting ansatz. By calculating and 
comparing their divergences, curls, and time derivatives, we improve the ansatz. Thus improved 
ansatz is further improved through the same procedure and the final ansatz fields are identified 
with Jefimenko's fields. Our inference of Maxwell's equations is in much the same spirit as the 
derivation of static differential equations (divergences and curls) from the Coulomb's and Biot- 
Savart's fields known experimentally without knowing their governing laws of the divergence and 
curl equations. 
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I. INTRODUCTION 



There have been incessant interests in the derivation of Maxwell's equations in the 
literature. ^"^^ The special theory of relativity, by providing the transformation properties 
of the electromagnetic field vectors, permits the derivation of Maxwell's equations from the 
field equations satisfied by a static electric field. Recent works by Davis^^ and Heras^^ pro- 
vided especially nice treatments. Davis developed a generalized Helmholtz theorem for an 
arbitrary localized time- varying vector function and showed that the mathematical relations 
commonly referred to as Maxwell's equations can be derived from that theorem. Heras 
formulated an existence theorem that states that, given localized scalar and vector time- 
dependent sources satisfying the continuity equation, there exist two retarded fields that 
satisfy a set of four field equations. These two formahsms are complementary. 

Though Maxwell himself succeeded in establishing the system of Maxwell's equations, he 
did not obtain the electric and magnetic fields in terms of the sources in general motion. 
Jefimenko's equations, ^^'^^ which are expressions for the electric and magnetic fields pro- 
duced by time-dependent charge and current densities, have received much attention only 
recently^^"^° and some books^^"^^ contain a section on Jefimenko's equations in their recent 
editions. 

The goal of this paper is to infer Jefimenko's and Maxwell's equations at the same time 
by adopting an appealing idea of retardation, which dictates the true nature of the gen- 
eration and propagation of the fields, as a basic axiom. The essential idea is to build up 
the time-dependent fields by starting from Coulomb's and Biot-Savart's fields, and system- 
atically improving them. Our inference of Jefimenko's equations, we believe, can be more 
accessible for the students of an intermediate electromagnetism course than derivations of 
those equations given in the literature. 

The organization of the paper is as follows: In Sec. II, we review briefly the continuity 
equations for sources with retarded time argument and summarize static fields and their 
differential equations. The starting ansatz with retarded integrals for dynamic fields is 
introduced Sec. III. The detailed process of improving previously taken ansatz is described 
in Sec. IV, and the identifications of Jefimenko's and Maxwell's equations are made in Sec. V. 
The last section is devoted to the conclusion. In the Appendix, we list some formulas of 
vector analysis, especially, involving retarded source density functions. 
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II. CHARGE CONSERVATION AT A RETARDED TIME 



The notions of charge density p and current density J arise from one reahty, i.e., the very 
existence of electrically charged particles in motion: 

p{v,t) = ^g,(5(r-r,(i)), (la) 
J(r,t) = J]g.^5(r-r,(t)). (lb) 

The conservation of charge is expressed locally by the equation of continuity at time t' , 

V'-J(r', + ^^ = 0, (2) 

whose time derivative yields another continuity equation 

The del operation (V') in Eqs. (2) and (3) pertains to the first argument (spatial argument 
r'), not to the second argument (time argument t') of the operands J and dJ/dt'. We note 
that Eqs. (2) and (3) are satisfied at any time t'. Thus, at a (field) retarded time r at a 
space point r' defined relative to (present) time t at space point r by 

T = t-\r-r'\/c, (4) 

the above continuity equations should take the following forms 

V'-[J(r',r)].fi.ed + p(r',r) = 0, (5a) 
V'-[j(r',T)],fi,ed + p(r',T) = 0, (5b) 

where the dot over p(r', r) and J(r',r) denotes derivative with respect to the second argu- 
ment, T. Note that df{r',T)/dr — df{v',T)/dt for any retarded field f{r',r), because |r — r'| 
in the definition of t is a constant independent of the time t. This remarkable point will be 
exploited later frequently. Because of the bizarre dependence of J(r',T) and j(r',T) on r', it 
follows that 

V • J(r', r) = -p(r', r) + - j(r', r) • , (6a) 

c |r — r| 

V • j(r', r) = -p(r', r) + -J{r', r) • , (6b) 
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by virtue of Eqs. (5a) and (5b). 

In passing, the notion of time-varying density fields given by Eqs. (la) and (lb) enables 
us to avoid a (particle) retarded time Tj endowed to each dynamic particle at the position 

rj(Ti) 

Ti = t-\r-ri{Ti)\/c. (7) 

As a matter of fact, due to this explicit and implicit dependence of Tj on the expressions of 
the Lienard-Wiechert fields for individual charges have very complicated forms. However, by 
simply superposing the Lienard-Wiechert fields due to the individual charges and introducing 
densities given by Eqs. (la) and (lb), one might reach Jefimenko's equations. 

In what follows, we will consider a system of charges in otherwise empty space; the 
discussion omits consideration of matter in bulk. In statics (p = and J = 0), it is well 
known that the Coulomb's electrostatic field produced by its source p, 

47reo Jv* |r-r'|3 
and the Biot-Savart's magnetostatic field produced by its source J, 



47reoc^ Jv* |r - r'l'^ 



satisfy 



V • E(r) = -p(r) , (10a) 

Co 

V-B(r) = 0, (10b) 

V X E(r) = , (10c) 

V X B(r) = ^J(r) . (lOd) 

In Eqs. (8) and (9), the integration region V* is all space (i.e., the entire universe). The 
electric and magnetic fields given by Eqs. (8) and (9) are characterized by the superposition 
principle through the spacial integrations. On one hand, we may regard the governing laws 

for fields [Eqs. (lOa), (10b), (10c), and (lOd)] as the consequences of direct applications of 
divergence and curl operations on empirically established physical entities [Eqs. (8) and (9)]. 
We may regard, on the other hand, the fields of Eqs. (8) and (9) as the solutions to the the 
governing laws. 



4 



In transition to the realm of dynamical fields, authors of conventional textbooks usually 
introduce, in historical sequence, the empirical law of Faraday induction first and then the 
theoretically deduced law of Maxwell induction. In the next section, we will attempt to 
modify directly the field expressions themselves with the idea of retardation as a guiding 
principle, and then the governing laws, if any 



III. A STARTING ANSATZ SET WITH RETARDED FIELDS 

Instantaneous action-at-a-distance extensions of Eqs. (8) and (9), i.e., fields Eo(r, and 
Bo(r, t) obtained simply by replacing the argument r' of p and J in the integrands of Eqs. (8) 
and (9) with (r', t) do not satisfy the continuity equation, differently from the static fields of 
Eqs. (8) and (9). We will not try other types of instantaneous action-at-a-distance extensions 
any more. 

Let us apply this time the idea of retardation to the fields of Eqs. (8) and (9) by replacing 
the argument r' of source functions with (r', t) to get the following starting ansatz set: 

To check whether these fields respect the law of charge conservation, we first single out p and 
J by taking divergence and curl on Eqs. (11a) and (lib), respectively. From the expression 
for the divergence of Ei, 



I j» j»/ 3 j» j»/ 1 3 



dV , (12) 



and Eqs. (A. 3) and (A. 6), we obtain the charge density 

p(r, t) = 60 V • E, (r, t) + ^ j^^ -^T^.^V . (13) 

We note that time derivative of this charge density, in general, does not vanish for time- 
varying sources, i.e., 

p(r,i)^0. (14) 
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In a similar manner, from the curl of Bi, 



V X Mi{r,t) 



V* 



r — r 



r — r' 



r — r 



_.^-V)J(r',r)-(J(r',r).V)- 



and Eqs. (A.3), (A.8), (A. 12), and (A. 15), we obtain 

J(r,t) = eoc^V X Bi(r,t) , 



(15) 



(16) 



where on the right-hand side a vanishing surface integral over the surface S* surrounding 
the volume V* 



1 f J(r', T){xj - x'j) 



r — r 



'13 



■dA' 



(17) 



47r 

has been omitted. 

Because the divergence of J of Eq. (16) vanishes identically, the continuity equation 
demand that the time derivative of p of Eq. (13) also should vanish identically, which con- 
tradicts Eq. (14), our original assumption of time- varying source. 



IV. IMPROVEMENT OF THE STARTING ANSATZ 

The origin of the dilemma arisen in the last paragraph of Sec. Ill lies in that the right 
hand side of Eq. (13) becomes p(r, t), which is identical to the left-hand side of the equation. 
Therefore, in order to get out of the dilemma, we should attempt to modify the starting 
ansatz electric field, first, by adding E2(r,t) such that 

which absorbs the last term on the right-hand side of Eq. (13) and allows an additional term 
r(r,i). Then Eq. (13) becomes 

p(r, t) = eo (V • [Ei(r, t) + E^ir, t)] - r(r, t)) . (19) 

Reminding two formulas in vector analysis, Eqs. (A. 2) and (A. 18), one can find 

^ ^ / '^':,:''':/^ dV' (20) 



47reo Jv* c|r — r' 
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to solve Eq. (18). By virtue of Eq. (A. 7), one can also find F to be fixed as 

p(r'. t) 



4.7reoJv* c2|r-r'| 



c^|r — r 

which can be expressed, by virtue of Eqs. (6a) and (6b), as 



r{r,t) 



V* 



1 •• V — r' 
V'-J(r',r)--J.^ 

c r — r 



1 



r 



-dV, 



(21) 



(22) 



and, further by virtue of an identity V'(l/|r — r'|) — — V(l/|r — r'|) and Eq. (A. 9) in reverse 
order, as 

1 r j(r',T) 



where we have omitted a vanishing surface term 

1 f j(r',r) 



dA'. 



Eqs. (18) and (23) give the divergence of modified electric field Ei +E2: 



dV 



(23) 



(24) 



(25) 



Since we are dealing with dynamic fields, it is useful to take the time derivative, in 
addition to curl and divergence, of the dynamic fields. In these calculations, in order to 
relate, if possible, the electric field to the magnetic field, we will express the results in terms 
of J and J by use of the charge conservation laws [Eqs. (6a) and (6b)]. By virtue of the 
interchangeabihty of two kinds of differentiations d/dr and d/dt upon retarded fields, and 
Eqs. (6a), (6b), (A. 16), (A. 17), (A. 8), and (A. 9) read in reverse order, we see that the time 
derivative of Ei + E2, 



|[E.(M)+E,(r,,)l = ji-^^^ 



p(r', r)(r — r') p(r', r)(r — r') 



r-'|3 



c r 



dV , (26) 



becomes 
d 



-[Ei(r,t) + E2(r 



V 



r — r 



./|3 



V-J(r',r) 



r — r 



./|3 



+ 



(j(r',r) 



V 



r — r 



c|r — r'|2 \ ' / c|r — r' 

On the right-hand side of Eq. (27), we have omitted a vanishing surface term 



V- j(r',r) 



r — r 



y|2 



dV. (27) 



1 



p/|3 



+ 



c r 



»/|2 



dA'. 



(28) 



Using Eqs. (A.3), (A.4), (A.12), (A.13), and (A.19), we eventually obtain 



4( 1 / |£>),V')-ij(M). (29) 

dtyAneo Jv* c^\r -r'\ ) eo 



Beginning with a correction term E2, we have ended up with Eq. (29), which is a natural 
modification of Eq. (16). It can be readily shown that thus modified equations for p and J 
[Eqs. (25) and (29)] do not violate the charge conservation law. 

Next, let us calculate other type of spatial variation (i.e., the curl) of the modified field 
El +E2. Relations between the spatial and time variations of the electric and magnetic fields 
which do not confiict with the continuity cquation(s) will be the governing differential laws 
what we are seeking. Using Eqs. (A.4)-(A.7), we see that the expression for curl of Ei + E2, 

V X [Ei(r, t) + E2(r, t)\ = ^ j^^ \^p{r', r) V x - x Vp(r', r) 

results in 

V X [Ei(r,t)+E2(r,t)] = 0. (31) 

The results up to now for the divergence, curl, and time derivative of Ei + E2 are sum- 
marized, respectively, in Eqs. (25), (31), and (29). 

V. IDENTIFICATION OF JEFIMENKO'S AND MAXWELL'S EQUATIONS 

Guided by Eqs. (25) and (29), we define E3 and B2 as follows 



1 f j(r^r) X (r -rQ 
Aneoc^ Jv* c|r - r' 



^2(r,^) - 7^ ^ WI2 ' dV. (32b) 



Then, Eqs. (25) and (29) can be expressed in much neater forms: 

V • [Ei(r, t) + E2(r, t) + E^ir, t)] = -p(r, t) , (33a) 

Co 

V X [Bi(r, t) + B2(r, t)] = ^J(r, t) + ^^[Ei(r, t) + E2(r, t) + E^{r, t)] . (33b) 
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To obtain an equation for the curl of the newly defined quantity Ei + E2 + E3, it should 
suffice to calculate V x E3 because we already got Eq. (31) in hand. Also, the calculation 
of the divergence of the newly defined quantity Bi + B2 will be carried out in due course. 
By virtue of Eqs. (A.l) and (A. 11), the curl of E3, 



V xE3(r,t) 



4ne 



Jv* 



Jfr',r) X V 



c^lr — r'l c^lr — r 



-Vx J(r',r) 



dV , (34) 



becomes 



V X E3(r, t) 



V* 



j(r', r) X (r — r') j(r', r) x (r — r' 



r — r 



'13 



+ 



c r — r 



/|2 



(35) 



which enables us to write Eq. (31) as 



d 



V X [Ei(r, t) + E2(r, t) + Es{r, t)] = -^[Bi(r, t) + B2(r, t)] . 



(36) 



By virtue of Eqs. (A.4), (A.5), (A.IO), and (A.ll), the divergence of Bi +B2, 



V-[Bi(r,0+B2(r,t)] 



r — r 

|r 



r — r 



-— ■VxJ(r',r)-J(r',r)-V^ 



+ 



r — r 



c r — r 



-^■Vx J(r',T)-J(r',r)-Vx 



r r 



c r — r 



/|2 



dV, 



(37) 



becomes 



V- [Bi(r,t) +B2(r,t)] = 0. 



(38) 



Main results of this section are summarized by four equations: Eqs. (33a), (33b), (36), 
and (38). It is well known that two source equations [Eqs. (33a) and (33b)] respect the law of 
charge conservation. We remark that E3 and B2 of Eqs. (32a) and (32b) have been dictated 
by El, Bi, and E2 of Eqs. (11a), (lib), and (20); E3 and B2 are contained in Eqs. (25) and 
(29). We emphasize that E3 and B2 should not be considered to have come from the already 
known answer, the Jefimenko's equations. 

Identifications of Ei + E2 + E3 and Bi + B2 with the physical electric and magnetic fields 
E and B, i.e., Jefimenko's equations, and of Eqs. (33a), (33b), (36), and (38) with Maxwell's 
equations complete our simultaneous inference of Jefimenko's and Maxwell's equations. Of 
course, the validity of these identifications should be confirmed ultimately by experiments. 
At the classical level, numerous experiments say that electrodynamics of Maxwell's equations 
(or equally Jefimenko's equations) equipped with Lorentz force law works very well. 
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VI. CONCLUDING REMARKS 



We have demonstrated, in detail, the process of inferring Jefimenko's and Maxwell's 
equations, starting from a starting ansatz set (retarded versions of Coulomb's and Biot- 
Savart's fields) for time- varying fields. The method presented here is based on the knowledge 
of (i) the electrostatic and magnetostatic fields and the continuity equation, with which 
most students are familiar at the intermediate level, and (ii) the seemingly^^ complicated 
del operations on the retarded fields. 

Of course, the inference of Maxwell's equations themselves here is more comphcated than 
the remedy by Maxwell himself of the pathology of charge non-conserving source equations. 
Nonetheless, the inferring process here has its own merit. It enables us to get Jefimenko's 
equations with ease comparatively. Indeed, to obtain solutions to Maxwell's coupled equa- 
tions except in the simplest situations is a very complicated task. For sources in arbitrary 
motion, the derivation of Jefimenko's equations is known to be quite complicated; the knowl- 
edge of the generalized Helmholtz theorem (or the wave field theorem called by Jefimenko 
himself) or the Green function technique is prerequisite. The inferring process here, how- 
ever, does not require either knowledge. In this respect, it is hoped that our inference will 
provide a useful tool in exposing advanced students as well as intermediate-level students 
to the Jefimenko's equations at early stage in the electricity and magnetism course. 



Appendix: Useful formulas 

We hst some vector-analysis formulas (i) not involving retarded fields 

1 r - r' 

r - r' 1 

^•1^3^= ^T-^' (A-2) 



r — r' 



V-^— ^ = 47r<5^(r-rO, (A.3) 
r — r' 

Vx^— ^=0, (A.4) 

Vxr^^^O, (A.5) 
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and (ii) involving retarded source density functions 



Vp(r',r) = -ip(r',r)^^, (A.6) 



c l'^ ~ , 

Vp(r',T) = -^p(r',r)^^, (A.7) 

V-J(r',r) = -ij(r',r)-^, (A.8) 

V-j(r',r) = -ij(r',r).^^, (A.9) 

VxJ(r',r) = ^j(r',r)x^^, (A.IO) 

Vx j(r',r) = ij(r',r)x^:i^, (A.ll) 
c |r — r I 

r^-v)j(r',r) = -i<^. (A.12) 
|r_r/|d j c|r — rp 

Jl^.v)j(r',r) = -to. (A.13) 
|r — / c|r — r| 

The formulas for del operations on retarded quantities in this Appendix may be unfamiliar, 

but they can be derived readily from a basic differentiation rule, 

df{r',T) df{r',T) dr 1 df{r' ,t) xt - x[ 



dxi dr dxi c dr |r — r| 

We also list other useful formulas: 



(A.14) 



-X- (J(r^r)j^) =^:f^V^J(r^r)-(J(r^r).V^ (A.16) 

(j(r^r)^LZ_l^ =^^V^J(r^^ (A.17) 

V-(/G)=/V-G+(V/)-G, (A.18) 

V X (F X G) = (V ■ G)F - (V ■ F)G + (G ■ V)F - (F ■ V)G , (A.19) 
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